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INTRODUCTION 


Shortly after publication of the article referred to, we had some 
doubts about the validity of the energy method used to determine the 
rate of damping in space from the rate of time damping. 

A more complete analysis of the problem seemed to show, in fact, 
that only first order terms (in v l ! z ) were appropriately calculated, 
whereas terms of second order (in v) contained possible errors. 

To finally settle this question and, at the same time, give the exact 
expression for the spacial damping coefficient, it was necessary to 
repeat entirely the calculation made for spacial damping. This quite 
long and delicate calculation, carried out by Carry, shows definitely 
that the second order terms given by the energy method were incorrect. 
Nevertheless, the error introduced in this way was negligible for almost 
all practical cases. For this reason La Houille Blanche , busy with other 
things, asked us to suspend a publication which is nonessential and 
repeats a previous one. However, the importance of the question was 
renewed with the publication by Miche of a work on the "Properties of 
Ocean and Laboratory Wave Trains" which, furthermore, was recently 
analyzed by La Houille Blanche . In a note in his work (page 93), Miche 
claims without formal proof that the principle term in spacial damping 
is in v and not in v l U 4 and suggests that the fact of having obtained a 
term in v 1 ^ results from the use of the energy balance method. 

The calculations of Carry which are independent of the energy 
balance method show that that is not the case and that our conclusions 
with respect to the existence and importance of the term in v l h remain 
unchanged; when the speed of the wave on the bottom becomes appre- 
ciable, this term predominates. The term becomes negligible only 
when the depth is greater than one-half the wave length and when, as a 
result, the speed on the bottom is less than a few percent of the speed 
on the surface. 

We believe that Miche 1 s error results from his omitting to write 
that the speeds go to zero on the bottom. This condition, as is known, 
corresponds to the usual hypotheses of theories of viscous fluid flow. 


1 



While Biesel has calculated the time damping of a wave, it is some- 
times useful to have the expression for spacial damping; the calculational 
scheme is the same, but we shall repeat it for the sake of clarity. 

We know that there exists a current function of the form i]j = \\)\ 4- 
where 1^1 and ^2 satisfy, respectively, 

Ai|i 1= 0 


and 


The excess of pressure over the hydrostatic pressure is given by 
the total differential, 


d p = p 


St 


Hi 


dz - 


Hi 


d x 


or with being the conjugate function of cpi > 
5 <S>! 


P = - P 


St 


+ C' . 


(1) 


For the case G>f waves propagating in a finite depth H, vjj can be 
written in the form, 

= (A cosh az + B sinh az + C cosh (3z 

, ^ • u a \ _iax + ibt 

+ D sinh (3z) e , 

where the x axis coincides with the free surface, and the z axis is 
perpendicular to it., (3 is determined from the relation, 

p 2 =o 2 H • < 2 > 


The aim of our study is tb find the complex numbers a and S, thus 
defining the function i|j. 

The conditions which must be satisfied on the bottom are 
Sib 

= 0 
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and 


5i)j 


= 0 


for 

z = - H. 

If, for simplification, we let 
L = cosh aH, M = sinh aH, 

P - cosh pH, and Q = sinh pH, 
these conditions can be written as 

(a) AL - BM + CP - DQ = 0 

(b) (AM- BL) a + (CQ- DP) p - 0. 

If r|(x) designates the ordinate of the free surface, one also has 

iax + ibt 


5 ~n _ /H\ 

6t \ 6 x/ z = 


o 


ia (A + C) e 


from which one finds 


a . . iax + ibt 

r| + — (A + C) e = 0 . 


Furthermore, the expression for the pressure is 

iax + ibt 


that 


p = - pgz - pb ( A sinh az + B cosh az) e 
At the free surface (z = r|), the normal pressure should be zero so 
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p+ 2pv = 0 
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(c) Ag + B (2a v + ib) + Cg -rg- -2Daf3v = 0. 


The condition of zero tangential force is written as 

S 2 ^ . ill _ n 

6x 2 6 z 2 


or, using Equation (2), 


(d) 2 Aq 2 v + (2a 2 v + ib) C = 0. 

Eliminating A, B, C, and D from the preceding equations gives 


2a 2 v 

L 

Ma 

go- 

ib 

0 

-M 

•\La 

2a 2 v + ib 

2a 2 v + ib 

P 

QP 

_ga_ 

ib 

0 

-Q 

-P(3 

2apv 


Solving this determinant, one gets 

-4qV(3 2 QM +(3 |pMga + PL (2a 2 v + ib) 2 + 4a 2 v 2 PL 

- 4a 2 v (2a z v + ib)] -ga 2 QL - QMa (2a 2 v + ib) 2 = 0. 

Dividing by PL and replacing L, M, P, and Q by their values, 


that 


- 4a 3 v 2 (3 2 tanh aH tanh pH + (3 j ga tanh aH + (2a 2 v + i b) 2 } 


+ 4a 4 v 2 - 


4a 2 v 


cosh aH cosh (34 


( 2a 2 v + ib)] - ga 2 tanh (3H 


one finds 


(3) 


- a tanh (3H tanh aH (2a 2 v + ib) 2 = 0. 



Equations (2) and (3) can be used to find a and (3; (3 2 , on the order 
of 1 /v, will be very large so that we can set tanh pH = 1 and cosh 
PH = oo. Thus Equation (3) becomes 

- 4a 3 v 2 p 2 tanh aH + p (ga tanh qH - b 2 + 4ib a 2 v 
+ 8 a 4 v 2 ) - ga 2 - a tanh aH (2a 2 v + ib ) 2 = 0. 

Solving Equation (2) for p by a series expansion, we have 

a _ n /^~ (1 + j) a 2 (1 - V 1 ^ 

P = y/2 vi/z + 4/b- 

Replacing (3 in Equation (4) by this expression, we get 

- 4a 5 v 2 tanh aH -i 4a 3 v b tanhaH -ga 2 - a tanh aH (2a 2 v + ib) 1 

- ^ Vl/2 j (tanh aH ga - b 2 + 4iba 2 v 

+ 8a 4 v 2 ) = 0. 



Neglecting terms of order greater than v l /z , this becomes 
(tanh aH ga - b 2 ) - ga 2 4 - a tanh aH b 2 


( 6 ) 


+ v 1 ^ 


4iba 2 \ZTT ^ + (tanh aH ga - b 2 ) 


= 0 . 


A solution for a is looked for of the form 


a — Ci + C 2 + C 3 v + . . . 

Noting that 

tanh aH = tanh CiH + — + 

cosm CiH 


C 3 H 


cosh 2 CiH 


and equating the different powers of v in Equation ( 6 ), 
equations: 


C 2 H 2 tanh CiH 
cosh 2 CjH J l 

we have the three 
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gCi tanh Ci ~b 2 = 0, 


-gCi + b 2 Cl tanh CiH + 


\/Ell + 0 


gC 2 tanh CiH + 


gCiC 2 H 
cosh 2 CiH 


and 


(i ' ' 2 % ClCz + b2 ' ( C 2 tanh c i H + — Cl ^ 2 ^ ' ) 

Vb~ (1 + i) f ^ (C 3 H - C 2 H 2 tanh CiH) C|H 

ST 1 8 1 cosh> C,H + 8 ToshS.cTiT 

+ gC 3 tanh CiH 1 = 0. 

As is known, the wave number a = 2 tt/B of waves of angular fre- 
quency b in a perfect fluid is given by ga tanh ah = b 2 . By analogy we 
shall use here the notation Ci = a, but it must be emphasized that a 
does not represent the same concept for waves in viscous fluids. 


Thus, we have 

c = (1 - i) a 2 

2 .\/b~(sinh 2aH + 2aH) 


and 


4ia 3 

b (2aH + 2 sinh 2aH) 

2 sinh aH cosh 3 aH 
(sinh 2aH 4- 2aH) 2 


sinh 2aH 4- 


1 

sinh 2aH 4- 2aH 


from which 


a = a 


2a 2 v 


(1 -i) 


(si 


x (sinh 2aH 4- 


sinh 2aH + 2aH 
1 


i4 


a 2 v 


b (2aH 4- 2 sinh 2aH) 


2 sinh aH cosh 3 aH 


sinh 2aH + 2aH (sinh 2aH + 2aH) 2 


>] 


(7) 


and 


=/ET 

\h 


(1 + i) + 


a 2 (1 


i) 



( 8 ) 



From Equation (7), the formula giving the damping law as a function 
of distance traveled can be written as 


2h = 2h 0 exp 


sinh 2aH + 2aH 



sinh 2aH 


+ 


1 


sinh 2aH + 2aH 


Q 

2 sinh aH cosh aH 
(sinh 2aH + 2aH)2 


! )l 


Before comparing this formula to that obtained by Biesel using the 
energy method, a preliminary remark must be made on the definition 
of a and b. 


In the case of damping in time and not in space, treated !, in extenso" 
by Biesel, the wave length, L, was defined unambiguously, allowing 
the determination of a from the formula, 


a 


2t t 

T7 


b was then defined from a: 
b = y/ ag tanh aH 

For the case of spacial damping, the concept of wave length is 
no longer obvious since the profile of the wave at a given time is no 
longer periodic in space, the amplitude decreasing from one crest to 
another, a could be defined as the real part of a, but we have preferred 
to proceed more rationally starting from the periodicity in time. The 
latter is free from ambiguity. Thus, we have taken 



and have defined a as the solution of 
b = v/ag tanh aH 

In this way, we avoid talking of a wave length that does not correspond 
to an obvious physical reality, as Biesel avoided speaking of a period 
in his article. 

Thus, the definitions of a and b cannot automatically be transposed 
from the time damping case to the spacial damping case; strictly, Biesel 
would have had to have given the definitions of these quantities at the 
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same time as the last formula of his article. Furthermore, note that 
the various definitions that can be logically considered differ by 
quantities of the order of v l h. As a result, to the order of approxima- 
tion used, only terms in v in the damping formula should be affected 
by the choice of definitions of a and b. 

Thus, the differences that can be noted between the terms in v in 
our formula and that of Biesel can be partially explained. However, 
there exists another possible cause for the difference in the estimation 
of the rate of energy propagation; the latter may be effected by the 
viscosity, the most important term in the expression for this effect 
being probably of the order of v 1 ^. Strictly, it would be necessary to 
understand this modification to perfect Biesel’ s energy method and to 
obtain correctly the terms of the order of v. 

Ultimately, if one wanted to obtain the terms in v using the energy 
method, it would be necessary to complicate the method and, in doing 
so, take away the simplicity that is its primary advantage. 

The comparison between our formula and that of Biesel is clarified 
by the preceding remarks. Such a comparison shows that the term in 
v 1 /* is correctly given by the energy method whereas the term in v is 
in error. 


The practical importance of the differences between the two for- 
mulas is, nevertheless, not great, since: 


1) For large relative depths the difference disappears. 

2) For moderate depths (less than a half wave length) the differ- 
ence in the term in v is small and, furthermore, the term v l U 
predominates. 

3) For very small relative depths the term in v becomes the most 
important, and the asymptotic values are different, the damp- 
ing being 


2h = 2h 0 exp 
instead of 

2h = 2h 0 exp 


8 (aH)‘ 




ax 


= 2h 0 exp 


3v 


8bH 2 3 


ax 


1 

(£A J 


V 

4 (aH)2 

l b ) ax j 

— ^ iIq exp 

4bH2 ax 


In any case, it is important to note that the method used cannot 
give sure results for these orders of magnitude, since the assumptions 
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made during the calculation imply that v is very small and that h is 
not too small (hypothesis P/Q = 1). But the numerical examples given 
by Biesel show that that difference is significant only for liquids with 
* viscosities much greater than that of water or for very small relative 

depths (on the order of 1/60 of a wave length). 

Thus, the general conclusions of Biesel’s article remain tenable, 
at least qualitatively. In particular, it is true that in depths that are 
not too large, the damping is linked particularly to the term in v 1 ^. 

This is in contradiction to the statement of Miche on page 93 of his book, 
Proprietes des Trains d'Ondes Oceaniques et de Laboratoire (Properties 
of Ocean and Laboratory Wave Trains), edited by COEC. 

To finish the calculation we shall look for the expression for the 
current function, i|j. Equation (d) gives 



Equations (a) and (b) give 

^ __ Aa + C [a cosh aH cosh (3H -(3 sinh pH sinh aH] 
a cosh aH sinh pH - p sinh pH sinh aH 


and 


B - Cp + a [P cosh pH cosh aH - a sinh aH sinh pH] 
p cosh pH sinh aH - a cosh aH sinh pH 

Inserting A, B, and D in Equation (1), we get 
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cosh pH 


P 2 / L 


cosh pH 


cosh PH 


Z 


sinh a (H + Z ) 



cosh aH sinh P(H+ Z) 
cosh pH 



tanh pH cosh a (H + Z) + 


sinh pZ 
Z cosh pH 


with 


1 / P \ 3 c a cosh pH 

1 Z \ a / p sinh aH cosh pH - a cosh aH sinh pH 

a and p being defined by Equations (7) and v 8). 


It should be noted that the form of the current function is the same 
for time damping if one sets a = a and ib = K, a and K being defined 
by Biesel, and where 


K 


ib { l ~ ' sinhW Vb ) 


+ b 


1 

sinh ZaH 


+ 2 


a 2 v 


cosh 4aH + cosh ZaH - 1 
cosh 4aH - 1 



and 





b_ 

2v 


(1 + i) 


'ZSlV 


2 sinh 2aH 


+ 


a 2 v 


8 b sinh 2 2aH 


x (1-4 cosh <2aH - 2 cosh 4aH) ( 


1 - i) . 
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